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The decay r — » nnniy is analysed using different methods to account for the resonance structure, 
which is usually ascribed to the ai. One scenario is based on the recently developed techniques to 
generate axial-vector resonances dynamically, whereas in a second calculation the Oi is introduced 
as an explicit resonance. We investigate the influence of different assumptions on the result. In 
the molecule scenario the spectral function is described surprisingly well by adjusting only one free 
parameter. This result can be systematically improved by adding higher order corrections to the 
iterated Weinberg- Tomozawa interaction. Treating the ai as an explicit resonance on the other 
00 ' hand leads to peculiar properties. 

o ■ 
o : 

CN . I. INTRODUCTION 

c : 

^ ^ ' The constituent quark model [J d, [1] has been very successful in describing part of the observed hadron spectrum, 
i especially for the heavy-quark systems, e.g. charmonia and bottomonia [1]. On the other hand, especially in the 
' light-quark sector, there is still a lively debate about the nature of many hadronic states. One sector with a lot of 
activity is, for example, the light scalar meson sector (cr, ao(980), /o(980), k(900)). These states can not be explained 
within the naive constituent quark model, and many models have been proposed to explain the phenomenology of 
Qh^ these resonances. The suggestions for the nature of these resonances vary between qq states, multiquark states, KK 
Q^, bound states and superpositions of them (see e.g. P, [E Hi and references therein). A different route to explain 
^ the low- lying scalars has been taken in [3, Q (see also references therein). In these works the authors explain the 
, states as being dynamically generated by the interactions of the pseudoscalar mesons. The scatte ring amplitudes are 
calculated by iterating the lowest-order amplitudes of chiral perturbation theory (CHPT) m Uil, which leads to 
I ' a unitarisation of the amplitudes and creates poles which can be associated with the scalars. 
^ ' A similar question about the nature of hadronic resonances one encounters in the baryon sector, where the quark 
' model also has trouble to describe the baryon excitations and their properties in a satisfying way (see e.g. [13, fl3| 
T-H , and references therein). As in the scalar case, an alternative approach to explain the resonance structure has been to 
generate resonances by iterating the leading order interactions of a chiral effective theory. The pioneering work in that 
'"-^ , direction has been done in [IJ, [15| and was followed by many other works [S [Sill, ill [2^, 

which suggest a number 



^""^ ■ of — i baryon resonances to be generated dynamically by the interactions of Goldstone bosons and baryons, e.g . 
A(1405) and iV*(1535). Studying the interaction of the pseudoscalar mesons with the decuplet of baryons [111, [22| 
also led to the generation of many known — | resonances, as e.g. the A(1520). 
■ Recent works applied the approach to the interactions of the octet of Goldstone bosons with the nonet of vector mesons 
.5^ ' focusing on the = 1+ sector [2^(23]. The authors calculate the scattering amplitude by solving a Bethe-Salpeter 
. equation with a kernel fixed by the lowest-order interaction of a chiral expansion. The leading-order expression for 
' the scattering of Goldstone bosons off vector mesons in a chiral framework is given by the Weinberg- Tomozawa (WT) 
term 25, 26] and leads to a parameter free interaction. The only free parameter in the calculation enters through the 
regularisation of the loop integral in the Bethe-Salpeter equation. Poles have been found, which have been attributed 
to the axial-vector mesons. 

A comparison of the pole position and width is necessarily indirect and depends on the model, which is used to 
extract these quantities from the actual observables. In addition, the height of the scattering amplitude, or in other 
words the strength of the interaction, is not tested in this way. In the following we apply the method of dynamical 
generation directly to a physical process, namely the t decay. The r decay offers a clean probe to study the hadronic 
interactions since the weak interaction part is well understood and can be cleanly separated from the hadronic part, 
which we are interested in. The r decay into three pions is dominated by a resonance structure, which is usually 
ascribed to the ai (see [1] and references therein). Many of the references in [l| are based upon a parametrisation 
in terms of Breit-Wigner functions^which leads to model dependent results. The relation of our calculation to some 
more microscopic descriptions [13, HI] is discussed below. 

We calculate the r decay in two different ways: We first calculate it by assuming that the ai is generated dynamically 
and use the method from [l^ll^] to describe the decay ('molecule scenario'). This means that in this framework the 
r decay is essentially described as follows: From the weak interactions a pair of mesons emerges (one pseudoscalar 
meson, one vector meson). Their final state interaction produces the resonant ai structure. This process is depicted 
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FIG. 1: (a) Basic diagram describing the dynamically generated a\ in the r decay and (b) additional diagram, when the a\ is 
included explicitly. (j> and V are the intermediate pseudoscalar and vector meson, respectively, which can be either up or KK* . 



in Fig. [T]Ja), where the blob stands for the iterated loop diagrams. There are at most two free parameters in that 
calculation (in the simplest scenario only one), which enter in the renormalisation of the divergent loop integrals. 
All other parameters are fixed by chiral symmetry breaking and the properties of the p. In a second calculation, 
we introduce the ai explicitly. Here the idea is that the ai is a quark-antiquark state. At the hadronic level this 
substructure is not resolved and the ai should be included as an elementary field. A similar approach u sing chiral 
effective field theory including elementary vector mesons and axial- vector mesons has been performed in [27|. This 
method yields a good description of the spectral function for the decay into three pions. However, the width of the 
ai in [27| has been parametrised, whereas we generate the width by the oi decay into Goldstone bosons and vector 
mesons. In [2^ the authors successfully describe the spectral function for the decay t~ 2'k'^tt~v in the framework of 
the linear cr-model. The width of the ai in this model is generated from the elementary decays of the ai. Nonetheless, 
there is still a fundamental difference between the approach we suggest in the present paper and the works 
Also for our second calculation with an elementary oi we still include the WT term since there is no reason to neglect 
it. The essential additional diagram is shown in Fig. [TJb), where the blob again represents the iterated loop diagrams, 
but that time the kernel also includes the ai interaction, which is discussed in detail in Section IIVI Having both 
calculations at hand - one with and one without an elementary oi - we can compare both to experiment and see which 
scenario is favoured by the data. Since there exist excellent data for the r decay [1^ , one can expect that the results 
will be quite decisive. In case that the first scenario is favoured by experiment, this would be a sign that the oi is a 
dynamically generated resonance (molecule state) and in case the second calculation is favoured, this would be a hint 
that the oi is a quark-antiquark state. Some of the results of these calculations have already been shown in 30]. In 
the present paper we present much more details of the calculations and show additional results, as for example the 
investigation of Dalitz plot projection data from [3l| within the molecule scenario. 

The oi is especially interesting, since it is considered to be the chiral partner of the p [O, [s^]. One expects a chiral 
partner for every particle from chiral symmetry. Due to the spontaneous symmetry breaking, one does not find 
degenerate one-particle states with the right quantum numbers. Nevertheless, the chiral partners have to exist, not 
necessarily as one-particle states, but at least as multi-particle states. Unmasking the ai as a bound state of a vector 
meson with a Goldstone boson would therefore approve its role of the chiral partner and disapprove its existence as 
a one-particle state. In the meson-meson and meson-baryon scattering examples, mentioned before, one can also see 
that some of the dynamically generated resonances would qualify as the chiral partners of the scattered particles, 
although the question of the chiral partner for these particles is not as clear as for the ai and the p. Even for the 
chiral partner of the p a different suggestion besides the ai exists, namely the 5i(1235) [33j . 

The work is structured as follows: We first discuss the general framework of our calculations in Section[ITl and we write 
down the relevant interaction terms, which we will need during the calculation. Next we describe the unitarisation 
procedure to account for the final state correlations in Section ITTTl In Section ITVl we calculate the matrix elements for 
the T decay in the different scenarios. Afterwards in Section |V] we compare our results to experiment and in Section 
IVII we give a summary and an outlook. Further details on the formalism can be found in the appendix. 



II. CHIRAL INTERACTIONS AT TREE LEVEL 



The low-energy dynamics of the Goldstone bosons is described by chiral perturbation theory (CHPT) [ol. fiol. [lll| . 
To lowest order the Lagrangian reads 
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with 

[/^ei0/-Fo^ ^^2Bo{s + ip) (2) 



where 



and the covariant derivative 



V2Tr- -TfO + -^r, \/2A'" | (3) 
V2K- V2Tf -^?7^ 



D^,U = df,U - tr^U + iUl^ . (4) 




r^, lfj_, s and p are external fields, which promote the global SU(3)xSU(3) symmetry to a local one. The interaction 
of the weak gauge boson with the Goldstone bosons can be determined by setting (see e.g. jTl|) 

r^=0, l, = -^iW+T+ + h.c.), (5) 
where h.c. refers to the hermitian conjugate and 

/O Vud Vus\ 

(6) 

The scalar field s incorporates the explicit chiral symmetry breaking through the quark mass matrix. In the following 
it is sufficient to use s = diag(m„, m^, mg), Fq is the pion decay constant in the chiral limit and Bq parametrises 
the connection between the quark masses and the Goldstone boson masses. Following [IJ we use Fq = 90MeV 
throughout this work. 

The Lagrangian EqjT]) contains the direct coupling of the weak current to three pions, which describe the r decay at 
very low energies [34 1. At higher energies, however, the r decay is dominated by resonances, most notably the vector 
mesons. In the molecule scenario we assume that the process is driven by the decay into Goldstone boson and vector 
meson. The structure, which is usually attributed to the ai, is generated by the strong final state interactions of the 
Goldstone bosons and the vector mesons. This means that we need the interactions of the vector mesons with the 
Goldstone bosons and with the W boson. 



A. Vector meson couplings 



Several models have been proposed to introduce the vector mesons in the chiral Lagrangian, e.g. the Hidden 
Symmetry approach [3^ or the WGGWZ [H, [s^,!!!] scheme. Most of them were motivated by the phenomenological 
successful ideas of vector-meson dominance and universal coupling. We will use the WGGWZ scheme, where these 
features are implemented by putting constraints on the couplings. Besides the choice of the scheme to introduce 
the vector mesons, one also needs to choose the interpolating fields for the vector mesons. Instead of descri bing the 
particles in terms of four- vectors, the vector mesons can also be represented by antisymmetric tensor fields |39l. |40| . 
The approaches are of course equivalent since the choice of fields can not influence the physics. However, due to 
the truncation in momentum the two descriptions can differ by higher order contact terms (see e.g. [4C), 41]), which 
especially influence the behaviour at higher energies. In this work we use the WGGWZ scheme and describe the vector 
mesons by vector fields. In addition, we include contact terms to improve the high-energy behaviour. 
The octet of vector mesons is given by 



V2p; -pI+uII^ I . (7) 



is an admixture of the physical states tj^ and (j)^ (for details see e.g. We do not care about the details of this 
mixing, since these states do not contribute to our calculation. 

One can define a convenient representation by introducing the auxiliary quantity u, which is the square root of U [Tlj 

= U . (8) 
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The transformation on U under the chiral group induces a transformation on u, which is given by 



u — >u' = VmJlJ = RUK-\L,R,U). (9) 

The transformation of the vector fields under the chiral group in terms of K is given by 

V^^V;, = KiL,R,U)V^K^iL,R,U). (10) 

The SU(3) matrix K carries the SU(3)^ x SU(3)j^ transformation in a non-linear way. The covariant derivative is 
defined for any object X, which transforms as 

\7^X = d^X + [T^,X], r^^^{u^{d,,-tr^)u + u{d^-il^)u^)). (11) 

In the following we also need V^i/, the field strength tensor of the vector mesons, which is given by 

V^,, = V^V, - V.F,, . (12) 

Using vector fields in the above representation the already mentioned Weinberg- Tomozawa (WT) term is contained 
in the kinetic part of the Lagrangian. Together with the remaining relevant couplings of the vector mesons to lowest 
order i one is led to the following interaction terms 

^vec = -^Tr[[y^ d,V.]n - :^Tr[y^./7] - !^Tr[VK, , (13) 

with 

Ui_, = i[u^di_, - irf^)u - u{d^ - il^)u^], 

and 

/^"^ = uF^^v) ± u^F^^'u , 

where -F^J^ are the field strength tensors of the external left- and right-handed vector fields 

F^'' = a^r-9''/^-i[;^,r], F^" ^^^'r'' -^''r^' -i[r'',r'']. (14) 

The first term in Eq. (|13p is the WT term, which is parameter free. In a heavy- vector formalism [23l . |4^ it is the 
only term of order 0{q^), where q is the momentum of the Goldstone bosons in a chiral counting. The Lagrangian in 
Eq. has already been written down in . We note that the definition of in differs from our definition by 
a factor of \/2, which yields different coefficients in front of our terms. The two parameters fv,gv can be determined 
from the decay of the p into dileptons and two pions, respectively [40], which yields 

0.154 GeV 0.069 GeV , , 

In [13] the authors also give a theoretical estimate for these parameters, which yields 

, V2Fo 0.127 GeV Fq 0.064 GeV , „, 

These values slightly differ from the experimental values, and in Section|V]we will study the influence of this difference. 
It will turn out that the experimentally measured values Ea. ([T5|) describe the data best. 

Transforming the vector-field Lagrangian Eq. (jl3p into a Lagrangian employing tensor fields, one finds that one can 
account for the difference resulting from the choice of fields by adding the following term (see e.g. [iol. l4ll. |43|. |43| ) 

Cho - -2Tr[j^,j'^n (17) 

with 

fv 

Thus, the entire Lagrangian we use is 

£ = £2 + ^vec + ^ho ■ (18) 

We will study the importance of Cho below in Section |Vl 
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FIG. 2: Iteration of loop diagrams, corresponding to the approximation to the Bethe-Salpeter equation by using the WT term 
as kernel. Full lines denote vector mesons, dashed lines Goldstone bosons. 



B. Axial- vector meson couplings 

In the second scenario, where we introduce the ai exphcitly, we also need to add the coupling of the axial-vector 
mesons to the Lagrangian. The nonet of axial- vector mesons is given by [24] 

/a; + /i(1285) V2a+ V2K+^ \ 

A^^i V2ar -a; + /i(1285) V2K^^ • (19) 
V V2K-^ V2K% 72/1(1420)/^ 

The additional Lagrangian we are going to use to describe the interactions of the ai is 

Ca.M = -^TiiA^^f^"] + iciTT[V''''[A^,u,]] + ic2TT[A'^''[V^,u,]] , (20) 

where the first term incorporates the coupling of the ai to the W and the last two terms describe the decay of the 
ai into Goldstone boson and vector meson with the unknown constants ci and C2. The first term has again already 
been written down in ^4Q |. whereas one can find different approaches in the literature in order to describe the ai 
decay vertex. In [45| the authors propose a phenomenological Lagrangian in terms of tensor fields, which is successful 
in reproducing the decay branching ratios. In (46 ^ the hidden symmetry formalism was used to derive the pertinent 
terms, which yields the same results as the phenomenological approach. Comparing the vertex resulting from the 
Lagrangian above to these works, we find agreement by choosing 

Ci = , C2 = — . 21) 

4 ' 8 ^ ' 

Below we will also study variations of ci and C2 around these values. Looking at Eq. ((20)l we see that both terms 
describing the decay into vector and pseudoscalar meson contain one u^. This means that integrating out the axial- 
vector fields would generate to lowest order an interaction term of vector mesons and Goldstone bosons, which contains 
two and therefore leads to an expression of order 0{q'^). Since the WT term is 0{q^) including both interactions 
is not double counting. This will be an important aspect for our results presented below. 



III. UNITARISATION PROCEDURE 



Unitarisation methods have been used to extend th e ap plicability of the chiral Lagrangians to higher energies, e.g. 
the N/D method [1], the inverse amplitude method [47| or partial summations using the Bethe-Salpeter equation 
[l^llj]. In 0] it was shown that the N/D method is equivalent to a summation of diagrams using the Bethe-Salpeter 
equation with the kernel taken onshell. We follow the work in [H, and use the Bethe-Salpeter equation to describe 
the final state interactions between vector and pseudoscalar mesons. In the molecule scenario this serves to generate 
the axial-vector mesons dynamically. In the scenario with the explicit ai the width of the latter is generated. The 
iteration (of a point interaction) is shown diagrammatically in Fig. [21 In case we are not including the ai explicitly, we 
will use the amplitudes from jlj, H^], i.e. the iterated WT interaction, in order to describe the final state correlations 
of the vector meson and the Goldstone boson. In the scenario, where we explicitly take into account the ai, we 
include, in addition to the WT term, the ai interaction (cf. Fig. [3]). This in principle generates the width of the oi 
by the decay into vector meson and Goldstone boson with the WT term as correction. 

Concerning the choice for the kernel in the case that we include the oi explicitly, we note the following. The kernel 
of the Bethe-Salpeter equation can consist of parts, which are analytic in the energy region, one is interested in, and 
parts which are non-analytic. The non-analytic parts are possible s-channel resonances, which have a pole in the 
physical region. The analytic part consists of contact terms, as well as t- and u-channel processes, since they do not 
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FIG. 3: Kernel of the Bethe-Salpeter equation when including the a\ explicitly; i.e. one has to replace the point interaction in 
Fig. [5] by the two diagrams on the right hand side. 



have any singularities in the physical energy region. The analytic parts can be expanded in powers of the momenta 
of the involved particles (and in powers of the Goldstone boson masses). This leads to contact terms, which can 
equivalently and systematically be expressed in terms of a chiral Lagrangian. The approach we take is to consider 
the relevant s-channel resonances in the kernel and keep contact terms up to a specific order. This approximation is 
a model assumption, and it is not guaranteed that it works for the quite large energy region, we are interested in. 
However, it is certainly worthwhile to study its properties. Following this strategy, one has to avoid double counting 
between the s-channel processes and the contact terms, since for s < m^es these s-channel resonance terms can also 
contribute to the analytic part. Restricting the contact terms to the WT term, there is no problem in our case, 
since the WT term contributes at 0{q) and the elementary ai at 0{q^). Thus, there can be no double counting. 
In particular, this means that only considering an ai and neglecting the WT term is very questionable. The focus 
in discussing the results will therefore be on the interplay between the ai and the WT term. In order to check the 
systematics of our model, we will also study the influence of keeping contact terms up to order 0{q^) instead of 0{q) 
in the scenario without explicit oi. 

We will briefly summarise the formalism, which we employ to iterate the respective kernels, which basically follows 
the same lines as in [2^. The Bethe-Salpeter equation for the coupled channel problem is 

T;t{q.q,w)=Kf^{q,q^w) + J2 J ■^^Kfp{qJ,w)Gf^{l,w)Tt{l.q,w) (22) 

where T is the scattering amplitude, K the kernel, q (q) is the incoming (outgoing) momentum of the Goldstone 
boson, w the total four momentum and 



Id 



cd * ^"'^ 

^ {w- ly -ml + teP- M2 +ie^"' ^^^^ 

is the two-particle propagator. The indices a, b, c, d indicate the channel, which in our case can be either irp or KK* . 
The indices are chosen such that denotes the scattering of (&, i^) {a,fi). 

We will solve the Bethe-Salpeter equation by employing an expansion in helicity amplitudes as follows 

K;:i<ir,,w)^ J: v;(-ir^(.)r-('^^'^^)(9,g,^), (24) 

J,M,P,i,j 



T'Jiqrq^w)^ Y: Mii^'^\s)Y:;-'^''''^\q,q,w), (25) 

where s — is the invariant mass of the hadronic final state. The indices i,j correspond to the helicities and can 
be either or 1 in our calculation. The objects y.^'^('^*^^) are similar to the projectors in [2^ and are discussed 
in Appendix [X] In order to arrive at an expansion with coefficients which only depend on s, one has to take the 
amplitudes onshell. This has been discussed in several works and we refer to 0, [231 . [45| for a justification. 
In [23| the authors were concerned that the projectors are not analytic outside the centre-of-mass system (CMS) 
and therefore performed a transformation to covariant projectors, which mix different helicities. We will not do that 
transformation since we only work in the CMS and without this transformation it will be easier to connect to states 
with definite angular momentum (see Appendix [B]) . The projectors fulfil the following orthogonality relation 

/ (^^A',t^a(9,^«^)G"''a,«^)>^AfA/^'(^,9,«^) = '5A.A3<5pP'<5jy'Wr,X,(g,q,u;)(-V) , (26) 
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with the divergent loop integral 



hvis) = / ,o^^4.... n2 ^2^,. ,2 ..2 , ... • (27) 



The relation only holds up to additional tadpoles, which are dropped. Using this relation the Bethe-Salpeter equation 
turns into an algebraic equation for the expansion coefhcients 

cik Mcbkj{—I<f)v) ■ (28) 

c k 

We introduce the renormalised quantity J0y(s, /i) 

Jwi^^l^) = Uv{s) - l4>v{f^) , (29) 
which depends on the subtraction point /i. In order to render Eg. ipSl) finite we substitute 

Ic/yvis) j0y(s,Ml) , (30) 

which introduces the a priori unknown parameter /ii. It remains to determine the coefficients Vctik in order to calculate 
the scattering amplitude. The relevant formulas are given in Appendix El 

IV. r DECAY 
A. The decay width 

Since the weak decay vertex is common in all diagrams (cf. Fig. [5] and Fig. Illl below). we can separate it from the 
hadronic information by writing the invariant matrix element as 

tM = C{s)S, ^g^'' - -j^j , (31) 

where we used the following abbreviations 

= v{p^)j^{l - J5)u{Pt) (32) 

and 



Cis) = J^] . C . (33) 



Wu denotes the hadronic tensor, which we will calculate in detail below. Gp is the Fermi constant, which is connected 

4V2M2 



to the weak gauge coupling by Gf = J 2 ■ The decay width is given by 



2Mr 

where d04 is the four-body phase space. We define the hadronic tensor 



dT=^^\M\'dct>^, (34) 



W^,= / iy;iy.d<^3, (35) 

J 03 



which by Lorentz invariance must have the following structure 



W^''{s) = W,{s)[9^'' -]+W2{s)—^. (36) 

Plugging in Eq. ((3T|l into Ea. (p4)) and expressing the width in terms of Wi and W2 we get for the differential decay 
width for the process — > 2'k^t:~Vt 

■\C\\M^^^sf(w2~Wr(l + ^\] . (37) 



ds 2MrS ' I ^ ^ V \ M^J J ' 
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Wi and W2 can then be calculated from the following two equations 



SWiis) + W2is) ^ W*-Wd<j)3, (38) 

W2(s)^- [ wW*wWd<l)3- (39) 

The factor of ^ taking care of the identical particles in the final state has been introduced in Eq.([37]), and thus it 
does not appear in Eg. ([55)1 and Eq. (|5^ . The longitudinal part W2 has no visible effect on the calculation, which was 
expected since it is proportional to as can be seen below. Therefore, we will drop it from now on. 
The quantity, we compare most of our calculations with, is the spectral function ai(s) for the decay r — > 27r'^7r~^'r. 
The complete spectral function A{s) is defined by 

A(s) ^ -—^Uris) , (40) 
s 

with the hadronic vacuum polarisation 

n^i. = + Hl (41) 

and 

n^. - i f d^xe'^^{0\TA^{x)A,{0)^\0) , (42) 



where T is the time ordering symbol and A^^ is the charged axial current 

Af^ ^ujf^^^d. (43) 
In terms of Wi the spectral function for the decay t — > 2TT'~'TT^iyr is given by 

= (44) 



B. The general picture 

One can not expect chiral perturbation theory to describe the r decay in the whole energy region (see [s^ for 
pure CHPT calculations), since the energies, which are involved are beyond 1 GeV and the decay is dominated by 
resonance structures. Including vector mesons at tree level in the calculation will certainly improve the calculation, 
but still one can not expect to find a satisfying description of the data due to the strong correlations in the final state. 
In particular, the vector mesons at tree level can not produce an axial- vector resonance. In Fig. 2] we see the spectral 
function calculated in lowest-order CHPT and by including vector mesons in comparison to data. The lowest-order 
CHPT calculation using the Lagrangian Eq.([T|) (cf. Fig. [S^, O3) can only describe the lowest data points, which 
are far below the wp threshold. The onset of the rise in the region 0.5 — 0.7 GeV^ is described much better if one 
includes the tree-level vector-meson diagrams Fig. [5t, [5jl. Nonetheless, the main bump in the data at about 1.5 GeV^ 
is clearly out of reach. The philosophy in the present work is that the final state interaction is dominated by the 
coupled-channel dynamics of the np state. That point of view is suggested by the improvement when we include the 
p in the calculation and the height of the amplitude, which we can see in Fig. 3) The deviation from the data for 
higher energies can then be explained by the increasing importance of the rescattering diagrams, describing the final 
state interactions (Fig. [5^,[S]:). We neglect further correlations between the pions, e.g. the diagram in Fig. [51 which 
we expect to have only a minor influence, since the tree-level result for the three-pion process is much smaller than 
the tree- level result for the process including vector mesons (see Fig. [J]). 

When we include the ai we expect, of course, also the oi to have a major influence on the result, but we recall that 
we still include the WT term in the calculation. 

In Section[V|we will also look at Dalitz plot projections, where one can see the clearly dominating p in the final state, 
which gives another justification for the choice of processes, which we include. 

In the present section we will present the formalism for different scenarios. Results are postponed to Section \V\ 
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FIG. 5: Relevant diagrams for the decay r~ 2-k n'v without including the a\. (j> and V correspond to intermediate 
Goldstone boson and vector meson {up or KK*). The blob represents the final state interactions obtained from the solution 
of the Bethe-Salpeter equation (see Fig. [Sjl. 



C. Calculation of r decay without a\ 



In Fig. [5] we see the processes which we take into account. The diagrams Fig. [S^, [SJd are the lowest-order CHPT 
processes, Fig. [SJ:, [5Ji are the tree-level processes including vector mesons and the diagrams Fig. [S^, [Sf describe the 
rescattering. from Eg. ([31]) is split into these contributions 



(45) 
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FIG. 6: Diagram describing pion correlations, which we do not include in our calculation. The blob denotes the final state 
interactions of the pions. 



where W^^ corresponds to the processes in Fig. [5^, Eb, W^^^. to the diagrams Fig. [5t, [5ji and Wjj^y to Fig. [5^, [5f. 
The first two functions are given by 



.rru ( gVudgV ( f ( 2 ( N N 

+ (/v - '^9v){{mi2 ■ qz)q2v ~ "2121.(93 ■ 92)) 



-((mi2 • g3)(w • 92) - (iw • mi2)(g3 • 92))^ ^72 — TT + (<?! ^ ^a) , (47) 



where qi and 93 are the momenta of the likewise non-charged pions, 92 is the momentum of the charged pion, 
w = 9i + 92 + 93, = qi + qj and the self energy of the vector mesons 11 is taken from [i^. {qi ^ q^) denotes 
the same amplitude with the pion momenta qi and 53 exchanged, which arise due to the appearance of two identical 
pions. The rescattering part is given by 

^ W^; + W^^, ^ (g^- - ^) br ^2^ ^l^, _ ^ (gi - 92). + (91-93), (48) 

where the nontrivial information is contained in bx- Wj^ contains an additional loop diagram, which needs to be 
renormalised. As in the Bethe-Salpeter equation ;23| we also drop additional tadpoles here, which leads to 



bT = ^-^^MnnJM ( ifv - 2gy)ks - ml + M^) + 2gy ( \mI + -^iK - K ^ 



AV2F^ ^^isii-^i^J^-lA^s) [{fv - 2gv)-{s - + M^.) + 2gv l^-Mf,. + ^^(m^ - M^. - s 



(49) 



The expansion coefficients for the scattering amplitude Matij have been given in Section IIIII We numbered the 
different isospin channels, where the channel 1 corresponds to irp and 2 denotes KK*. Thus, e.g. Mim denotes 
the amplitude for np scattering with the helicity of the p equal 1 in incoming and outgoing channel. We note that 
the renormalised loop integral J0y(/i2) appearing in 6t does not have to depend on the same subtraction constant 
as the loop integrals in the scattering amplitude. Therefore we denote the subtraction constant of this loop with p2 
and the subtraction constant appearing in the scattering amplitude with pi. We will discuss the appearance of two 
subtraction constants and their relation in more detail later. 

The scattering amplitude we used for this calculation was already determined in [2^ and [l^l and in terms of the 
expansion in the projectors from Appendix [XI we have 

Va\n=gab, (50) 

KbOl = /TT ° gab , (51) 
V 2Mva 
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FIG. 7: Higher-order contribution to the decay r 2n"n v. The blob with label '4' denotes contact interactions emerging 
from Eq. (|17|) . In a chiral counting these contact interactions are of fourth order. 



with 



9ab = [is - {Ml, + Ml, + + M^,) - -^{M^, - Ml,){Ml, - M^J ) , 



(54) 



(s + ~ M^J , 



(55) 



Pcma = 2:^.y/(s - [Mva + M^aY){s - {My a - M^a)^) 



(56) 



and 



WT 



2 -%/2~ 
-%/2 1 



(57) 



The scattering amplitude is then determined by Eq. (|^5|) . The second term in Eg. ([55)1 is in principle a higher-order 
term and we neglect it for the moment. We explicitly checked the influence of this term and there was no visible 
difference in the results by including the term. With these coefficients the lowest partial wave of the potential takes 
an easy form 



^ — QabYllfiU + 9ab 



V2Mva 



^1+ , ^b 
^Ol^v + 9ab 



V2Mvb 



. ^a^b 

" J 1 f 



(58) 



Choice of interpolating fields 

The matrix element above was calculated by using vector fields as interpolating fields for the vector mesons. In order 
to improve the high-energy behaviour we include additional higher-order interactions, which account for the difference 
in using vector fields or tensor fields [iO] . The additional terms are contained in Eq. (fTT)) , which leads to the following 
contribution 



iMho = CS^ ( g^"- - ) , 



(59) 



with 



w''w''\ / gVud9v 



ifv - 25y)(g2;y(m-i2 • qs) ~ mi2u{q2 ■ qs)) 



9Vud9vfv 



F§ s{s - ml) 



w^'{{w ■ q2) {17112 ■ qs) - {q2 ■ q3){w ■ mi2)) + {qi qa) ■ 



(60) 
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FIG. 8: Sum of diagrams, which build up the full p propagator. 



\ \ \ \ \ 

\ \ \ \ \ 

FIG. 9: Sum of higher-order diagrams, contributing to the t decay. 



which is shown diagrammatically in Fig. [T] We will show that the net effect of adding this contribution (in a resummed 
way) is the replacement 



n2 Ml 

^ ^ (61) 



p2 _ M2 - n p2 _ 2 _ n 



in W^ec fEa. (|T7)l V At a first look, it seems like we actually have to replace + 11. In order to show the 

validity of Ea. (j6ip . we look at the sum of diagrams building up the full p propagator, which are shown in Fig. [51 So 
far we have just included additionally the two diagrams, which can be seen in Fig. [71 But looking at the sum in Fig. 
[HI it would be reasonable to also include the sum of diagrams in Fig. [9l which we call Wi^tm- The central relation in 
this problem is that the higher-order contact term is proportional to the lowest-order diagram including the p, which 
is stated more precise in Fig. 1101 This relation guarantees that one can split off the pmr vertex in the higher-order 
contact term in the same way as for the resonance diagram. Thus, using the relation in Fig. IIOI both sums together 
yield 



=W^ \l - = W — - (62) 

* vec ' * * surn ^ * vec 1 2 1 vec 9 ' / 




which leads exactly to the replacement advocated in Eq. fCTl) . 

We omitted the diagrams with the tt intermediate state in the discussion (right diagram Fig. [7]), since the arguments 
follow exactly the same lines. 



D. Calculation of r decay with explicit a\ 



Next we include the ai explicitly. We introduce it as a bare resonance and generate the width by summing up the 
self-energy contributions from the decay of the ai in pseudoscalar and vector meson. This is automatically achieved 
by using the Bethe-Salpeter equation. As already mentioned, we also consider the WT term. Altogether, this leads 
to the diagrams shown in Fig. [11] In comparison to Fig. [51 there are two additional diagrams, where the W merges 
into the ai. Furthermore the blobs indicating the resummation are not the same as in the calculation before, since 
we add the ai interaction to the kernel. The additional process considered in the kernel is shown in Fig. [3l and it 




FIG. 10: Relation between higher-order contact term and lowest-order resonance diagram. 
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leads to the following expression 

^--^--^f73M|:l^"'^- MlT; (63) 
(ci((?>" - p ■ + C2(«/^g" - q ■ wg^'')){c^{r^ - p ■ qg^^^ + c^iw^q^ - q ■ wg^^^)) , 

where Mq^ is the mass of the ai, which is considered to be a free parameter. The expansion into projectors is given 

by 

Va\n = -^F,ab + gab (64) 

V^L = -^;^(^-^-« - ^-.x/ip^) + ^^Jab (65) 



•^0522 — 
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with 

■f^i = ^^ZY^ TfT ("^1 • 9) (P • 9) + C1C2 (g • w) (p • q) + C1C2 (p • g) (g • w) + C2 (g • w) (g • w) ) (68) 

F3 = ^^^r^ \.2 ((^1 - • g + (ciC2 - c^)g . u;) (69) 

Fi = ^^^r^ \7^{{cl-ciC2)p-q+{ciC2-cl)q-w)). (70) 

We recall that the coefficient Cwt is a matrix due to the coupled-channel structure of the problem. The indices a, 6, 
which were attached to the functions Fi in Ea. (|64ll70)) correspond to the respective channels. We also note that the 
Fi are the coefficients defined in Eq. l|A3p . The whole matrix element can be written similar to Eq. ([3T|l 



^M'^'=CS,[g^'^-^^W., (71) 

but this time with 

= W'JP + W'J^'^' + Wl" + W^^" + W'J"" + W^^'^P + Wf'^^' . (72) 

The last two terms correspond to the diagram Fig. [TT| :. W^^^'^ to the process in Fig. [TTh and W'J^'^ differ from W^"^ 
because of the different kernel. The sum of the modified contributions leads to 



(73) 



'12 ^"p 

The tensors L'^'^ and are given by 



5^ , —w^q ~ w^w . (74) 

s s 

The coefficients Ai, which incorporate the nontrivial part, are given in Appendix [C] The rest of the calculation follows 
the same lines as in the calculation before. We only have to substitute W^' with W^j^^^^ in the calculation of Wi. W2 
does not change, because neither WJ^ nor W^^^]^ contribute to W2 and we leave it out again. 

At the beginning of this section we used that the ai vertex, which results from the Lagrangian Eq. (|20p . is given by 

Ki = -^^^ci(0^ - P ■ qgn - ^^^C2{w'^q^ - ^ ■ qgn ■ (75) 
Using e^{p)p'^ 0, s = ^^a^i -^ai = 2Mp and ci = 2c2 the vertex can be cast into a different form 

Kl - ^4^^2(2gV -[2p-q + w q)g'^'^ + w'^ q>^) = ^^^C2(u; V " w ■ pg^ , (76) 

which is the actual vertex, which is used in |4^. We see that only with the simplifications above, the vertex of [i^ 
is the same as the one we use. These simplifications, however, basically mean to put certain momenta onshell and 
apply Weinberg's relation between the p and ai mass [s^ , as well as a relation between ci and C2 , which are anyway 
free parameters. We will discuss the influence of the difference between Ea. ([75|l and Ea. ([7S)) in Section IVBl 



E. Calculation of r decay including higher-order terms 



This time we again assume that the oi is generated dynamically. In addition to the WT term, we consider higher- 
order corrections to the kernel of the Bethe-Salpeter equation. One part of the higher-order correction, as the WT 
term itself, is contained in the kinetic part of the Lagrangian and it leads to the new contribution 



Kr ^^{q'-r-q^r)- (77) 
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Next we write down all terms with two pion momenta, which one can construct, taking care of parity, C-invariance, 
hcrmiticity and of course chiral symmetry 

+ A^iTriy^.y'-^u^T/^] + A^Tril^^.y^^u^u^] + y^Tx^^.u^^^'^u^ + + A[,Tr[y^u'']Tr[l/'^u,] 

+ A'ioTr[l/^u'^]Tr[Kii'^] + A'iiTr[V;u,]Tr[y''u'^] + A;,Tr[y^„u'^]Tr[Cu'^] + A'i3Tr[V;„u,]Tr[y''"u^] . 

(78) 

This together with Ea. ((77| leads to the following kernel in addition to the WT term 



P (^(A'l - 2A^)(g . q)g^^'^ + [X', - 2Xi, - + (A^ - 2X^ + \)rq'' - [K + - 2A^)(u; ■ q){w ■ q)g^^'^^ 

- ^((A[,(g • q) + iX[, + X[,){w ■ q){w ■ qW + X^^q^r + KxTq") ■ 

(79) 

We see that there are only eight independent variables contributing to the process. Dropping terms, which will not 
contribute to = 1+, we are down to six independent variables, which we call Ai, A2, . . . Ag. Thus, the kernel can 
be written as 

K: = ^(i^iil • 9) + A2(u; • g)(^« • q))9''- + Xsq^q") 

(«o) 



-^((A4(g -ID + X^iw- q){w ■ q))g^'' + AggV) 





Together with the WT term the expansion of the kernel reads 

Kihll = ~^Flab + gab (81) 

2 uj 

K621 = ~ o /77.^ {Flabi^a ~ Fiab^fspl) + /- ° Qab (82) 
KW2 = 3^jV,f^^ (-^b^lafc + F^abplV^) + -^^^^9ab (83) 



^622 = oAf ^Af (-^aUJbFlab - {XlCwTab - Xi5ab)^plpi + FsabUJbplVs + F4abUJaplVs] 



UJaiOb 

-gab- 



2MvaMvb^ 

(84) 



with 



-^1 = --T^ilo%)i^^ + «-^2) - ^('7ogo)(A4 + sAs) (85) 



p2 vi"-iuy\ ^ ■ pz 
^0 ^0 



_ Ch/t q-w 12 q-w 

-fs - - T^^— Ae (8b) 



_ CVt9- w I2 g-w, , . 

^A-—E^ ^A3-T^— — Ae- 

The diagrams, we have to include are the same as in Fig. Owith a different scattering amplitude describing the final 
state correlations, which leads to 



iM% = CS^ [ .g'^" - ^ ) , (88) 



with 



where 



VT^ = W'^'^P + W'^^^' + W^"" + Wf'' , (89) 



gVudgV 1 I \( Wf,Wa \ If 4,V , ON 0V\ Tl°' 

c^vJ^vKM \ gi,a j I \gvai + ^(/v - ^9v)al I L{ 



+ {gv4'' + \{Iv-'igv)<^ Lr 



(90) 



'12 



(qi - 92)7 + {qi ^ qs) , 
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/ / / 





FIG. 12: Diagrammatic form of the equation to determine the form factor of the W boson. The dashed hnes represent pions, 
the sohd hnes the p mesons and the wiggly hne the W boson. The bare vertex includes also the diagram with the intermediate 
pion (see Eg. ((93)) '). 



The coefficients can be found in Appendix O They are the same as in the case of including the exphcit a\ except 
that one has to replace the expansion coefficients of the scattering amplitude MaUj ■ 

F. W form factor 

Instead of first calculating the scattering amplitude, one could introduce the W form factor to determine the decay. 
Leaving out some details and only considering the WT term, it is possible to vifork out the decay width in a few lines. 
It is instructive to look at this simple calculation, since here the intermediate steps are not clouded by lengthy algebra 
and the core of the calculation is better visible. 
'Leaving out details' means 

• neglect longitudinal part of the hadronic tensor proportional to to^ 

• /y - 2gy = 

• neglect lowest-order CHPT diagrams (direct three-pion decays) 

• only 7rp channel, no coupled channels 

The first simplification actually has no visible influence on the result. The second simplification is numerically almost 
fulfilled, i.e. 

tLl^^Q.X ^ !v»{W-2gy) (91) 
iv 

The third approximation will only influence the very low-energy region of the decay. The only serious simplification 
is the last one, which we will later also see to have a minor influence on the results. Therefore, we can expect this 
slimmed down version to be pretty close to the full calculation. 
The equation, determining the W form factor, is 

V^^{q,w)^yr{^M+ J ^V'^"{l,w)G^^{l,w)K0''il,q,w) , (92) 
which can be seen in pictorial form in Fig. 1121 With the simplifications above, V^^'^ is given by 



where p is the momentum of the vector meson. We drop the term proportional to p'^ , since it will not contribute due 
to the form of the ptttt vertex (~ (gi — q2)^) and the renormalisation scheme, in which tadpoles are dropped. Thus, 
we get 

with Li;" defined in Eq. (jA4| . The kernel Rf^" is already known to be (see Eq. fSS]) ) 

if^'' = KoL^;;'' (95) 
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with 



4F2 



(96) 



Looking at Eg. ([92)1 and the form of the kernel, i.e. that it does not depend on q, we can write down a reasonable 
ansatz for V^^ 



^ V{s)L'{ 



Plugging in this ansatz in Eq. (f92|) we get 
and we can easily read off V to be 

V = 



12M2s 



{ml 



Ml 



^1 ' 



Vn 



1 



(97) 



(98) 



(99) 



The result is rendered finite by substituting /^rp —* J^pifJ-i)- The above calculation seems to employ only one 
subtraction point fii . This is in contrast to our derivation in Section IIV C[ where we argued that two different 
subtraction points can appear, namely one to renormalise the Bethe-Salpeter equation and one for the entrance loop 
from the W boson into the rescattering process. We will show in the following how the second subtraction point can 
also be recovered in the present calculation. 

Omitting the Lorentz structure for the moment, the full W decay vertex can be written as 



V = Vo + VoGT. 



(100) 



In the solution of the Bethe-Salpeter equation for the scattering matrix T, we will use G' in the following, in order to 
indicate a possibly different subtraction point. Using T = (1 — KG')~^K the form factor V can be written as 

V = Vn + VoG{l - KG'y^K = Vo{l - G'K){1 - G'K)-^ + V(,GK{1 - G'K)"^ 

= Va{l-G'K + GK){l-G'K)-\ (101) 

which corresponds to an equation of the form shown in Fig. I12[ provided one takes the bare W form factor as 

V^^Vo{l-G'K + GK). (102) 

We see that in the vertex we can effectively include a change in the subtraction point of the first loop relative to 
the subtraction point of the Bethe-Salpeter equation. We note that the change in the subtraction point is at least 
one order higher in a chiral counting, since the kernel K is already 0{q). We recall that Eq.((T3]) contains only the 
lowest-order W — > 4>V vertex. Using different renormalisation points for the loops G and G" gives us the possibility 
to account for modifications of this lowest-order expression. 

For the calculation of the whole decay, we only use one diagram, which is shown Fig. [T31 Thus, we get 



Wf 

* ' jorra 



with 



= 

jorrn 



iVg\ 



'-12 



M2 



n 



^1 (qia ~ q2a) + {qi ^ qs) ■ 



(103) 



(104) 



We do not show the results of the simplified calculations explicitly, since one can anticipate the outcome by looking 
at the discussions in Section |Vl In particular, in Fig. [20] we will see that neglecting the strangeness channel does not 
have a big effect. The other simplifications have already been estimated above to be less important. 
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TT 



FIG. 13: Diagram describing the r decay in the simplified version. 



0.25 




s[GeV^] 

FIG. 14: Spectral function for the decay r~ — > 2-k°'k~i't calculated with different choices for the interpolating fields, vl only 
uses the vector field Lagrangian Ea. (|13|l . whereas v2 additionally includes the contact terms from Ea. (|17|l . 



V. RESULTS 

A. Calculation without ai 

First we want to investigate the spectral function for the decay r~ 2n^TT~Vr calculated by iterating the WT 
term in order to dynamically generate the ai. We will discuss the influence of different aspects of the calculation on 
the results in detail and determine the values of the subtraction points. 

• Influence of interpolating fields and spectral distribution 

We discussed the different possibilities to describe vector particles, namely in terms of vector fields and in terms of 
antisymmetric tensor fields. We introduced higher-order corrections in order to account for the difference stemming 
from the choice of fields. For the present calculation we note that using the antisymmetric tensor fields, leads to the 
appearance of less derivatives, and therefore we expect a better high-energy behaviour. Instead of explicitly using the 
antisymmetric tensor fields, we use the vector representation but also include the higher-order contact terms given in 
Eq. lfTTl) . Fig. [Ml shows the spectral function calculated with vector fields (vl) and with vector fields including the 
higher-order terms (v2). One clearly sees the better high-energy behaviour for the case of v2. In Fig. [T3]we used 
dimensional regularisation with /ii = /i2 = Af^, which is the value from [2^. We will discuss the infiuence of the 
subtraction points in detail below. The kink which can be seen at about 1.9 GeV^ results from the threshold of the 
KK* channel. Using spectral distributions for the vector mesons, taken from [i^, smoothes the curve, which can be 
seen in Fig. [131 The curve also gets a little broader and moves to the right, but the overall structure is unchanged. 
If we do not state otherwise, the following calculations will always contain the higher-order corrections (i.e. v2) and 
the spectral function for the vector mesons in the loop. 
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FIG. 15: Spectral function for the decay — > 2-K''^iT^Ur calculated with and without including the width of the vector mesons 
in the loop integral. The curve labelled v2 is the same as in Fig. 1141 
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FIG. 16; Spectral function for the decay t~ — > 2n^'Tv~i/T calculated with different renormalisation schemes in comparison to 
data from [2911. 



• Influence of renormalisation 

In Section IIVI we encountered two subtraction points in our calculation. In the following we will investigate the 
influence of these two parameters. We start by setting /ii = ^2 and vary them simultaneously. In Fig. [11] we see the 
spectral function for the decay t~ — > 2Tr'^Tr~ Ut calculated with different renormalisation points in comparison to data 
from [29|]. The lowest curve is calculated using dimensional regularisation with = ^2 = Mp, which corresponds 
to the value employed in [2^. For the curve, which turned out to be the highest, we choose a cutoff scheme with a 
cutoff at 1 GeV in momentum, which corresponds to the choice in • The result in between was again calculated 
with dimensional regularisation, but this time with a subtraction point at 1^1=^2 = 8.5M^. Since including the 
spectral function for the vector mesons in the cutoff scheme is more complicated and does not lead to new insights 
for this comparison, we did not include the spectral distribution for the calculations shown in Fig. 1161 The different 
curves clearly differ in the height of the peak, but the position of the peak is not influenced. The width of the peak 
turns out to be too small in all prescriptions. 

It is also instructive to look directly at the scattering amplitudes, which describe the rescattering. We note that the 
scattering amplitude only depends on /xi and is independent of /X2. In Fig. [T7]we see the real and imaginary part of 
the scattering amplitude for irp scattering (corresponding to Mim) for different renormalisation descriptions. The 
curves shown in Fig. [T7] correspond to the highest and lowest curve in Fig. 1161 In addition, we plotted the lowest 



v2 

w.spec.fun. 




Hl = H2 = Mp"2 

= = 8Mp"2 

Cutoff 1 GeV 

data • 
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FIG. 17: Real (left) and imaginary (right) part of the scattering amplitude for yrp scattering. The curves correspond to the 
highest curve in Fig. (THland to the two curves shown in Fig. 1151 where the curve labelled 'v2' in Fig. [T^] corresponds to the 
curve labelled 'Dim. Reg.' in this picture. 
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FIG. 18: Spectral function for the decay r~ 2-K^-R~v.r, calculated by varying the subtraction point ^2 of the first loop and 
keeping the subtraction point in the scattering amplitude fixed at /ii = Mp, in comparison to data from [29j . 



curve using a spectral distribution for vector mesons in the loop, which corresponds to the result in Fig. [151 We see 
that the scattering amplitude hardly shows a resonance structure by using the subtraction point at (full curve in 
Fig. [T7)) . The resonant structure is more pronounced for the cutoff scheme (dashed curve in Fig. [T7)) . Including the 
spectral function of the p (dotted curve in Fig. I17p basically smoothes the curve. The bump in the imaginary part 
is moved to the left in the cutoff scheme, whereas in Fig. [11] one could hardly see a difference in the position of the 
peak. This shows that it is not so obvious to translate the structure seen in the scattering amplitude to the spectral 
function of the t decay. In other words, interferences between the tree level diagrams and the rescattering diagrams 
(cf. Fig. [5]) play an important role. 

In the following we will only use dimensional regularisation to render the loops finite. With this restriction we 
nevertheless cover the full discussion on the renormalisation parameter since there are only marginal differences by 
using different schemes, provided that the parameters are properly chosen. In principle this can be seen in Fig. I16[ 
where one can easily imagine that a further increase of the subtraction point leads to the same result as the one which 
was calculated with the cutoff. 

Next we want to investigate the effect of changing ^2 while keeping /ii fixed. We will use /ii = Af^, which in [2^ 
was determined by using crossing symmetry arguments. Thus, using these arguments to fix one subtraction point, 
we are in principle left with only one free parameter. The results for different choices of can be seen in Fig. 1181 
We see that we can describe the data very well by varying only that subtraction point and keeping /ii = fixed in 
the scattering amplitude [SO^. We note that choosing the subtraction point at ^2 = 8.5Mp, which describes the data 
best, corresponds approximately to a cutoff of 1 GeV in a cutoff scheme. Obviously this value is very reasonable. 
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FIG. 19: Spectral function for fii = Mp and fi2 ~ 8.5Mp (cf. Fig. llSp split up into the different contributions from the 
diagrams in Fig. [5] in comparison to data from [2^. The 'lowest order CHPT' curve corresponds to Fig. [5K,b, 'inc. vector 
mesons' to Fig. [5j;,d and 'only rescattering' to Fig. [5^,f. 



In Fig. [T9]we see the spectral function for /ii = Mp and /i2 = 8.5M^ (cf. Fig. [T8|) split into the different contributions 
from the diagrams shown in Fig. [51 We see that the bump is partly created by the negative interference of the 
rescattering diagrams and the diagrams including the vector mesons at tree level. The little bump, we see in the 
rescattering contribution alone appears at the wrong position and only the sum of all diagrams gives the pronounced 
peak, which is of course the only quantity that can be measured. 

• Influence of coupled channels 
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FIG. 20: Spectral function for the decay calculated with and without including the kaon channel in comparison 

to data from [29| . 

The spectral function calculated with and without including the strangeness channel is shown in Fig. 1201 which shows 
that the bump also appears without the kaon channel. The height is a little less by leaving out the kaons, but that 
could be compensated by varying the subtraction point /i2. The rise in the data in the energy region up to about 
1.1 GeV^ can also be described by leaving out the kaons, but the width of the peak is better described by including 
both channels. However, the effect is pretty small and one can safely say that the 7rp channel plays the dominant role. 

• Varying gv and fv 



data 
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lowest order CHPT 

inc. vector mesons 

only rescattering 
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coupled channels 

w/o kaon channel 




22 




s[GeV^] 

FIG. 21: Spectral function for the decay — + 27r'V^!/T calculated with different values for gv and fv in comparison to data 
from [2y|. 'gl' corresponds to the theoretically motivated values of fv and gv according to Eq. (|106P with the subtraction points 
/ii — Mp and fi2 = 8.5Mp. 'g2' uses the same parameters as 'gl', except fi2, which is chosen to be /12 ~ UMp. 



So far we used the experimentally measured values for fv and gv, which are given by 



0.154 GeV 0.069 GeV , , 

As already noted in Section |TT1 these values slightly differ from 

which are the values, obtained by theoretical considerations and approximations in [40| . In order to see the influence 
of these parameters on the results, we show in Fig. [21] the spectral function calculated with different values for 
gv and fv- 'gl' corresponds to the theoretically motivated values of fv and gv according to Eq. (|106p with the 
subtraction points /.ti = and /i2 = S.SAf^. 'g2' uses the same parameters as 'gl', except for /i2, which is chosen 
to be ^2 = 14M^ in order to fit approximately the height of the peak (and still get roughly the shape). We see 
that the moderate difference in fv and gv has a sizeable impact on the results, which is due to the fact that the 
combination fvgv appears quadratically in the final formulas. The difference of fvgv between using the values of 
Ea. (|105p and Eq. (|106p is about 25%. The change in the height of the peak can be compensated by a readjustment of 
the subtraction point /i2, but the spectral function in this case seems to be shifted to the right. We note that there is 
no other parameter, which potentially can influence the spectral function up to about s w 0.7 GcV^, as can be seen 
from the discussions before. 

Except of the influence at low energies and the resulting small shift, varying fv and gv seems to have a similar effect 
as varying /i2. This is not too surprising, since changing fv and gv changes also the W decay vertex and leaves the 
scattering amplitude untouched (cf. discussion in Section lIVFp . 

• Stable p 

In Section HVBI we discussed which diagrams we should include in our calculation, and we assumed that the contribu- 
tion from the pion final state interactions is small. In order to show that this is a reasonable assumption, we compare 
our previous calculations with one, where the p is assumed to be stable. This means we look at a spectral function 
obtained from the final state irp instead of Svr. With the notation from Section IIVI and neglecting the longitudinal 
part proportional to mj, we get 

^stable = ^[9 — j ^T~'^^^P> + [3 — J 2Fo 9.^o.~PaP^)ep{p) , (107) 
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FIG. 22: Spectral function for the decay — > -k"^ Vr calculated by assuming the p to be stable ('stable p') and spectral 
function for the usual three-pion final state ('pion final state') in comparison to data from [2^. In order to see the mere 
difference by assuming different final states, we do not use a spectral distribution for the vector mesons in the loop. The 
subtraction points are chosen according to the best choice at p\ = Mp and p2 = 8.5Mp. 



and therefore 

3 2(27r)5^V 9v 4i^o 9v ^ J ,,„„^ 

For simplicity, we also neglected terms ^ [fv — "^gv) in the tree level diagram {W^^^^). We checked explicitly, that 
these terms influence the results by less than 10%, and therefore they would only lengthen the formulas above. The 
negligible influence is of course expected from Eq. ipTj) . In order to see the net effect of assuming a stable p, we do not 
include the spectral function for the vector mesons in the loops. When the p is stable, the threshold for the decay of 
the T is of course sharper and moved to the right, which can be seen in Fig. [511 We see that beside these differences, 
the structure is the same as before and there is not much room for a big contribution from pion final state interactions. 

We want to summarise shortly what we saw so far. In a scenario where the ai is generated dynamically we 
employ two parameters pi and fi2- Varying both simultaneously with fj,i = H2, we saw that we always got a peak 
at the same position with varying height and a too small width. Using the value from [23j to fix p,i, we describe 
the data very well by choosing the only remaining free parameter p2 at S.SM^. We also saw that in the process 
under consideration the main contribution came from the irp channel, while the kaon channel plays a minor role. In 
addition, we investigated the infiuence of the parameters fv and gv and found that the results are quite sensitive to 
these parameters. In particular, the low-energy behaviour is best described by using the values, which are directly 
determined from experiment. Finally we assumed the p to be stable, which yields a qualitatively similar result. This 
eliminates concerns about possibly large final state interactions of the pions (cf. Fig. [S]). 



B. Calculation with explicit ai 



Now we want to look at the results of the calculation when we include the ai explicitly. A very small coupling 
of the ai together with the values from the calculations before will of course reproduce the results from before and 
will give a good description of the data. To check whether a scenario with an explicit ai can also describe the data 
reasonably well, we have to demand that the coupling is not almost zero, and we expect the value of the coupling to 
be comparable to the values found in [13, IH,!!^, i.e. comparable to Eq. ((2T|) . In order to get non zero couphngs and 
still be reasonably close to the data, we have to keep the contribution from the WT term small. We start by choosing 
fJ-i = fJ-2 = Mp according to [1^, which we can expect to be a good choice by looking at Fig. [121 Scanning through 
the parameter space, it turned out that in most cases a two bump structure is observed. We show an example of this 
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FIG. 23: Spectral function for the decay — > 2Tv°iT^Ur including the ai with different sets of paramaters in comparison to 
data from _29|. The parameter sets are given in Tab. ID 
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TABLE I: Different sets of parameters for the calculations with explicit ai. The remark 'w/o WT' means that the WT term 
is not included, the additional remark 'Eq. (|76p ' means that the ai decay vertex from Eq. (|76p is used (which does not employ 
the parameter ci) and 'fv,gv Ea. (|106p ' means that we choose the values from Eq. (|106p for fv and gv- 

in Fig. [22] (set 3). By finctuning the parameters, it is possible to merge these two bmxips into one, vifhich can also 
be seen in Fig. [521 (set 1 or set 2). The parameters leading to these curves are given in Tab. [H From set 1 and set 
2, we see that there are different possible choices for the parameters, which can (more or less) describe the data. A 
deviation is only seen for 0.8 GeV^ ^ s ^ 1-1 GeV^. 

Next we want to see, how far we get by switching off the WT term. We perform such an analysis without the 
WT term for completeness. We recall our strategy discussed in Section [TTTl to approximate the (j>V scattering kernel 
by possible resonances (here the oi) plus contact terms of lowest order. In that strategy there is no justification to 
neglect even the lowest-order contact term, which is just the WT term with its strength fixed model independently 
by chiral symmetry breaking. In Fig. [Mj we see the spectral function calculated with explicit oi but without the 
WT term in the kernel (set 4). In that case the second bump disappears and by changing the parameters, one can 
determine the position and the height of the peak. Although one might expect that the width of the peak can be 
adjusted by the choice of ci and C2, this is not the case, since there is a more complex interplay between the position 
and the width of the peak. There are lots of parameter choices which give a qualitatively similar curve. The curve 
labelled 'set 4' represents the best fit to the data by varying all parameters and we see that it agrees with the choice 
for Cl and C2 from 46] and 45] (cf. Eg. lj^T]) ). For the curve labelled 'set 5', we used the ai decay vertex from Eq. ljTS]) 
and we see that in this case the shape of the peak is described very well. This shows the uncertainties in the shape of 
the width and how it is influenced by the energy dependence of the ai decay vertex. Therefore, we have to be careful 
to draw conclusions from the exact shape of the width. 

From Fig. [Mj (set 4) and Fig. [221 one should not conclude that the WT term is a correction, which improves the shape 
of the width, since one does not simply switch on the WT term in order to get from Fig. [Mjto Fig. [221 Instead the 
parameter sets are quite different (cf. Tab. [l| and one has to finetune the parameters in order to arrive at a single 
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FIG. 24: Spectral function for the decay r~ 2ti n^Vr including the a\ with different sets of paramaters in comparison to 
data from [20|. The WT was not included in these calculations. In addition, the curve labelled 'set 5' uses a different energy 
dependence to describe the a\ decay (see text for details). 




reasonable peak, when one includes the WT term. 

Looking at 'set 4' in Fig. [21]one might worry that we do not have the most sophisticated model describing the explicit 
fli and that other models describe the data better (e.g. [13, H^). However, our description of the ai is completely 
sufficient to show the strong influence of the WT term. 

Next, we want to see the role of the WT term played in the best result, which was shown in Fig. [221 (set I). In Fig. 
[25] we plotted the result from parameter set I with and without including the WT term. Although with this choice 
of subtraction points the WT term is suppressed very strongly, it obviously still has a major influence on the result. 

In the discussion of the results without the explicit ai, we showed results for the theoretically motivated values 
of fv and gy given in Eq. (|106p . Using these theoretical values, it is possible to further suppress the WT term in 
comparison to the explicit ai. In Fig. [26|we see the best result, which we found in this case (set 6). In Fig. [211 we 
found that the rise in the data for energies between about 0.7 GeV^ ^ ^ < 1.2 GeV^ is described worse. Here, we find 
the same problem in the case of the ai. However, we note that using the theoretically motivated parameters for fv 
and gv it is much easier to get rid of the second bump. In this case one can hardly call it finetuning to obtain a single 
reasonable one-peak structure. Nonetheless, we stress again that the r decay data favour the experimental values for 
fv and gy given in Eq. (|105p as compared to the theoretically motivated values of Eq. (ll06p . Using the latter the rise 
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FIG. 26: Spectral function for the decay — » including the ai using different sets of paramaters in comparison to 

data from [291 1. 
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TABLE 11: Different sets of parameters, which yield a good description of the spectral function. 

of the data in the energy region 0.7 GeV^ <s< 1.2 GeV^ is underestimated (cf. Fig. [H] curves gl,g2 and Fig. 1261 set 
6). 

Looking at the results including the ai it is not so easy to draw an immediate conclusion. For sure, one can say 
that the WT term has a major influence on the result. The second bump structure can be recovered in almost every 
calculation including the oi together with the WT term. An important point is that the inclusion of the WT term 
leads to very strong effects, although we already kept the contribution very small. Only by finetuning one can merge 
the two appearing bumps. However, merging two bumps by finetuning the parameters does not seem to be a natural 
way of reproducing the data. In other words: Why should an elementary state appear right at the position, where an 
attractive potential already created a peak? We note again that the strength of the WT term is model independently 
fixed by chiral symmetry breaking. Since the WT term alone already produces a peak at the right position, one 
could already expect that a description of the data including the ai has to be accompanied by a delicate choice of 
the parameters. Still, it would be too much to talk about a definite sign that there is no explicit ai. However, the 
peculiarities with explicit ai together with the success of the description without the oi (molecule scenario) should 
be regarded as a good indication. In the next section, we will show that adding higher-order corrections to the WT 
term it is possible to systematically improve the situation in the molecule scenario and that the ordering of diagrams 
makes sense in this scenario without an explicit ai. 

C. Higher-order terms 

In Section lIV El we determined the corrections to the kernel at 0{q^), which led to six new unknown parameters. In 
the following we leave out the explicit oi again and show the infiuence of these corrections on the results. In Fig. [27| 
we show the spectral function with and without including the higher-order correction. There are several parameter 
sets (Tab. |TT]), which can describe the data in a qualitatively similar way. We see that the higher-order terms can 
be chosen such that they systematically improve the agreement with the data. Note that the size of the higher-order 
terms is not constrained by chiral symmetry (except that they should be of natural size - a demand of every effective 
field theory). 



27 





data * 






only WT 


- 




hoset 1 








- 






- 

























0.5 1 1.5 2 2.5 3 

s[GeV^] 



FIG. 27: Spectral function for the decay r 2ti -k Vr including higher-order terms in the kernel in comparison to data from 
[2^. For the choice of the parameters see Tab. Ull 



data 




s[GeV^] 



FIG. 28: Spectral function for the decay r~ 2n tt^Ut including higher-order terms in the kernel in comparison to data 
from 29]. For the choice of the parameters see Tab. [Ill In the calculation including only the WT term we employ the same 
subtraction point /i2 = 9 Mp as in the calculation including the higher-order terms. 



Next we investigate the connection between the higher-order terms and the subtraction points. Changing 112 can 
hardly be compensated by the higher-order terms, which is expected, since /12 acts as a higher-oder correction to the 
W decay vertex and not to the scattering amphtude. However, a slight raise in /12 can be compensated, as can be 
seen in Fig. [551 There we use fi2 = 9M^ for both calculations. We recall that without higher-order terms the best 
value was ^2 — 8.5AIp (cf. Fig. \TE\i . It can also be seen that the higher-order terms do not touch the energy region 
up to about s ~ 1 GeV^. 

Next we want to discuss the connection to /ii. Here we can expect that the higher-order terms are at least to some 
extent able to compensate for changes. In Fig. [29] we can see that we can account for moving the subtraction point 
/ii by changing the parameters of the higher-order terms. The parameters are again given in Tab. [HI We sec that the 
compensation is even better than expected, since the corrections up to order 0{q'^) do not carry all the structures, 
which might be influenced by moving the subtraction point. 

It is also interesting to look directly at the changes induced in the scattering amplitude. In Fig. [3D| we see the real 
and imaginary part of the scattering amplitude for the different parameter sets given in Tab. [TT] in comparison to 
the calculation without higher-order corrections. In order not to overload the figures, we show four different plots. 
Figs. [SOTa) and (b) show the scattering amplitudes for the first three parameter sets in comparison to a calculation 
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FIG. 29: Spectral function for the decay — > 2tv°tv^Ut with and without higher-order terms using ni = fi2 = 8.5Mp in 
comparison to data. The parameters for the higher-order corrections are given in Tab. [Ill 




FIG. 30: Real (left) and imaginary (right) part of the scattering amplitude for np scattering with and without higher-order 
corrections. The upper two plots show the scattering amplitudes for hoset 1-3 in comparison to a calculation without higher- 
order terms and /ii = Mp. The lower two plots show the scattering amplitude for hoset 4 in comparison to a calculation 
without higher-order corrections and /ii = Mp and fii — 8.5Mp, respectively. 
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without higher-order corrections and fii = M^. We see that the scattering amphtude is not modified much. Only 
'hoset 2' shows a different structure, which becomes most obvious in the imaginary part. This change in the scattering 
amphtude seems to be correlated with the parameter A3, which we will further investigate, when we look at the Dalitz 
plot projections in Section fV Dl Figs. [HOTc) and (d) show the scattering amplitudes for 'hoset 4', which was chosen 
to compensate for the change in the subtraction point /ii. The figure shows the scattering amplitude for /ii = 
and jii = 8.5Mp without higher-order terms in comparison to jii = 8.5Mp with higher-order corrections. We see 
that the corrections bring the scattering amplitudes for /ii — 8.5M^ back into the shape they had before, when we 
used Hi = Mp without higher-order corrections. In other words, changes in the renormalisation point can be replaced 
by changes in the higher-order terms. The renormalisation scale dependence is reduced as it should, when including 
higher-order terms. 

One might argue that including the higher-order corrections was unnecessary and describing the data with 6 
parameters is no success. However, the point in including the higher-order corrections is not that we can describe the 
data with seven parameters, but that they systematically improve the result. In case of the inclusion of an explicit ai, 
we saw that adding the WT term to the ai interaction worsened the results. Here, however, adding the higher-order 
terms to the kernel behaves as a correction. Note that the calculations are not ordered according to usual perturbation 
theory. Instead the kernel of the Bethe-Salpeter equation is calculated in perturbation theory. The convergence of 
that kind of perturbative expansion is not guaranteed. Therefore it is encouraging to see that the next to leading 
order terms behave as a correction and are even able to improve the agreement with the data. 

We note that there are many possible choices for the parameters, which describe the data. In Section IV Dl we will see 
that the different sets can be further discriminated by looking at the Dalitz plot projections. 

D. Dalitz plot projections 

In Fig. |3T]we show the Dalitz plot projections in 771^2 or rn^^ for the calculation using the WT term only and for 
the calculation including the higher-order terms in comparison to data from [sij. We determined the normalisation 
of the theoretical curve such that the area under all curves, corresponding to different slices of -y/s, agrees with the 
area under the data points. We also subtracted the contribution to the data, which was identified as background in 
[sH . Using this normalisation, we do not really lose any information, since we already saw before that the spectral 
function was well reproduced for all invariant masses, which implies a proper total decay width and therefore a proper 
normalisation. Fig. [31] clearly shows that the final state is dominated by the p meson, and the data are described 
quite well by all parameter sets. The last two plots seem to show an improvement by including the higher-order 
corrections. The improvement in the last two plots is most pronounced for 'hoset 2'. But since we overshoot at lower 
771^2 ^^'^ the error bars are pretty large for these invariant masses, the advantage is not very stringent. We recall that 
qi and ga are the momenta of the likewise pions and that the amplitude is symmetric under the exchange qi (73. 
Thus, 

771^2 s-ud 77123 s-re the invariant mass of the intermediate p, which we clearly see in the Dalitz plot projections. 
In Fig. 1221 we plot the number of decays versus 777^3, which is the invariant mass of the likewise pions. These pions 
do not build up the p and therefore the structure in these plots is completely different from Fig. |3TJ When we look 
at Fig. 1221 we also find that the calculations including the higher-order corrections describe the data better, which 
again is more pronounced for 'hoset 2'. The steep rise at small mfa is much better reproduced by 'hoset 2' and 
also the additional structure at higher invariant masses is reproduced better with 'hoset 2', although we overshoot 
that structure. We note that using 'hoset 3' and 'hoset 4', we do not get a noteworthy difference in the Dalitz plot 
projections in comparison to a calculation without higher-order corrections. Thus, the improvement in the Dalitz 
plots seems to be correlated with the parameter A3, which is non-vanishing for 'hoset 1' and 'hoset 2' (cf. Tab. HI)) . 
It is interesting to see the amount of d-wave contributions from the different parameter sets. In Fig. |33|we plotted 
the absolute value of the ratio of the respective coefficient Css and c^s to the sum of all coefficients (see Appendix 
IB]) . We see, that the parameter sets, which describe the Dalitz plot projections better, clearly have a higher d-wave 
contribution. Thus, our calculation indicates a population of d-waves in the r decay in the amount, which is shown 
in Fig. 1331 We note that the leading order contribution to Csd is given by the terms proportional to A3 and Ag in 
Eq. (|80p . A statement about pure d-wave transitions given by Cdd would be more complicated, since terms of higher 
chiral order than q^ would contribute at leading order. 

VI. SUMMARY AND OUTLOOK 

We calculated the process for different scenarios. In the first scenario (molecule scenario )we analysed 

the decay based on the recently developed techniques to generate axial vector resonances dynamically |23l . [23 | . The 
picture we promote is that the process is dominated by 7rp final state interactions, which are described by iterating 
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FIG. 31: Dalitz plot projections in or 77123 with and without higher-order corrections in comparison to data from fsH]. The 
different parameter sets can be found in Tab. |IT1 The curve labelled = 8.5Mp corresponds to a calculation using the WT 
term only, fii — Alp and fi2 = 8.5Mp. 
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FIG. 32: Dalitz plot projections in mfs with and without higher-order corrections in comparison to data from [sj]. The different 
parameter sets can be found in Tab. [Ill The curve labelled /X2 = 8.5Mp corresponds to a calculation using the WT term only, 
m = Mp and fi2 = 8.5Mp. 
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FIG. 33: Population of angular momentum states expressed in Css/{cs3 + 2csd + Cdd) (left) and Csd/{css + 'iCsd + c^d) (right). 
The curve labelled 'WT term' shows the population for the WT term only, whereas the other kernels include higher-order 
corrections with the parameters given in Tab. |TI] 



the WT term. The weak decay is part of the standard model and the WT term is predicted parameter free from 
chiral symmetry. The remaining couphng constants {fv,gv), which describe the interaction of the vector mesons 
are determined by the properties of the p (40| . The only unknown parameters in the calculation enter through the 
renormalisation of the loop integrals. We introduced two subtraction constants to render the loop integrals finite. One 
subtraction constant (^i) renormalises the loops in the scattering amplitude describing the final state interactions. 
This parameter was already introduced in and fixed by crossing symmetry arguments. The other subtraction 
constant {^2) enters in the renormalisation of the first loop, which contains the W decay vertex. We investigated the 
influence of these parameters on the results. First, we varied /ii and ^2 simultaneously with /ii — fi2 and compared 
in particular the different values used in [2^ and . We found that all choices produce a peak at the same position 
with a different height. The position of this peak was roughly in the region of the resonant structure seen in the data, 
but the width always turned out to be too small. Afterwards we investigated the influence of 112 by keeping /xi fixed. 
Using the crossing symmetry argument from '23'| in order to determine fj,i, leaves us with one free parameter. Fitting 
this parameter /i2, we reproduced the spectral function for the decay t~ — > 27r*'7r~z^,- quite well. 
In a second scenario, we explicitly introduced the ai in the calculation. This introduces new parameters, namely the 
mass of the oi, its coupling Ja to the W boson and the couplings ci and C2 to the vector- meson Goldstone boson 
states. The most obvious feature in that calculation is, that due to the strong influence of the WT term a second 
bump appears. Finetuning the parameters, one can merge the two bumps into one and the data can be described 
more or less satisfyingly. However, the results of these calculations are unsatisfying. An important point is that the 
inclusion of the WT term leads to very strong effects, although we already kept the contribution very small. Merging 
two bumps by finetuning the parameters does not seem to be a natural way of reproducing the data. Since the WT 
alone already produces a peak at the right position, one could expect already that a description of the data including 
the ai has to be accompanied by a delicate choice of the parameters. In addition, one can obviously not talk about a 
small correction, which is induced by the WT term. 

A further improvement of the molecule scenario was found by introducing higher order corrections to the kernel. 
These corrections introduce six new parameters and many combinations of these parameters could be found, which 
fit the spectral function very well. The correction induced by these terms were well behaved. That was not clear 
from the beginning, since we calculated the kernel of the Bethe-Salpeter equation perturbatively, which does not 
automatically guarantee that we picked up all important contributions for the scattering amplitude itself. Therefore, 
this is an encouraging fact, which puts further foundation to the calculation and shows its systematic nature. 
Comparing our calculation to the Dalitz projections, we found that including higher order terms, which carry d-wave 
components, describe the data better. We compared the size of the coefficient Cgd, which describes the transitions 
from an s-wave to ci-wave state, for different parameter sets. The size of this coefficient was clearly correlated with 
the qualitative description of the Dalitz plot data. 

To summarise, one finds that without the explicit ai one has a well behaved model, which can be systematically 
improved and which describes the data very well. Most parameters (in the simplest scenario, all but one) are fixed by 
chiral symmetry breaking and the well known properties of the p. Including an explicit ai leads to peculiar properties, 
if one tries to generate the width consistently from the Bethc-Salpeter equation and includes the WT term. When we 
tried to describe the data with an explicit ai the strength of the WT interaction caused the most severe problems. 
On the other hand, this strength is fixed by chiral symmetry breaking. In addition, we recall that taking into account 
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an explicit oi and the WT interaction is not double counting. Essentially we claim that the WT interaction should 
not be disregarded as has, however, been done in many previous approaches. On the other hand without an explicit 
fli the WT interaction has the right strength to generate a resonant structure dynamically. These indications point 
towards a dynamical nature of the ai as a (coupled-channel) meson molecule. 

As an outlook we note that a further step in the calculation would be to include medium effects in order to see 
what happens to the ai in case we approach the chiral symmetry restoration [ioj ]. In principle, when the restoration 
happens, the axial- vector spectral function, defined in Section FlVAl must be degenerate with the corresponding vector 
spectral function. In the latter the p meson prominently appears, at least in the vacuum [29j . It is, however, not so 
clear what chiral restoration implies for the specific part of the spectral function with a three-pion final state. In any 
case one would expect a drastic reshaping of both the vector and the axial-vector spectral function. 
It would also be interesting to figure out how well the molecule scenario agrees with QCD lattice calculations [sO] of 
the axial- vector current-current correlator (in the specific region accessible by lattice QCD). Here one has to perform 
the calculations with a higher pion mass in order to connect to lattice QCD calculations. This also brings into play 
pion mass corrections to the involved coupling constants as for example Fo,fv, gv [5l|. Still one can expect that the 
presented framework offers enough predictive power to obtain a valuable comparison to lattice QCD. 
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APPENDIX A: PROJECTORS 



We briefly summarise the most important formulas, which define the projectors and which are used to determine the 
expansion coefficients. A scalar amplitude for the scattering of vector mesons with a scalar particle can be expanded 
as follows [52} 

et (p, A)r'"'e,(p, A) = ^ —^Dj*^^i^, 0, -<P)Di,^{<f>, 6, -<P){J M\\T\J MX) , (AI) 

JM 

where p, A (p. A) are the momenta and helicities of the incoming (outgoing) vector mesons, D are the Wigner rotation 
functions and \ JM\) denotes a state with total angular momentum J, its projection M and with the helicity of the 
vector particle being A. Choosing the incoming particles to fly along the 2:-axis and the scattered particles to move in 
the xz plane, the formula reduces to 

el{p,\)T^'^e,(j>,\) = Y.^2J +l){X\T'\\)dl-^{e) , (A2) 
J 

where we omitted J, M in the denotation of the states and d are the Simplified Wigner functions or d-functions. By 
Lorentz invariance the scattering amplitude can be written in terms of five scalar functions Fi 

T^. = Y.F,V^,, (A3) 

i 

with 

' - X (A4) 

q — w 



where q (q) is the incoming (outgoing) momentum of the Goldstone boson. We note that there are only five independent 
terms, since terms containing p'^ or p'^ vanish due to e^{p)p^ = 0. Thus, using the orthogonality relation of the d- 
functions, one can express the expansion coefficients of Ea. (|A2p in terms of the Fi 

(A|T'^|A) = i^\t(p,A) (Y^F^lA e,ip,\)dl-^smede. (A5) 
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We further introduce parity eigenstates, which are given by 

(l±l = ^((-l|±(l|)- (A6) 

The defining equation for a projector with total angular momentum J, its projection M, parity P and helicities A2, A3 
is 

e^Hp,^i)YlXU^'l''yil^^^) - S^,,^,J,,^,^^^^^ (A7) 

For = 1+ and the kinematics described above the explicit form of the projectors is 

1+ _ 3 / 1 2 ^^2; 3 -ZU 4 



Y,^L.^M^[-—Ll,. + L 



2 \ ^ pps p\/s ^ p y's 

3 / ujx .2 , r4 1 



3 / LOX ,9 1 



(A8) 



SMMxLl^ 
pps 



where M, lu (M, ZU) are the mass and the energy of the incoming (outgoing) vector meson, pcm the centre-of-mass 
momentum and s = (p + q)"^ the total invariant energy of the process. For practical calculations in the centre-of-mass 
system, however, it is enough to know Eq. (jA7[) and the expansion coefficients from Eq. ljASp . 



APPENDIX B: CONNECTION BETWEEN HELICITY STATES AND ORBITAL ANGULAR 

MOMENTUM 

In order to determine the s- and d-wave component of the vector-meson Goldstone boson two-particle state, we 
need to know the relation between the helicity states and the orbital angular momentum I. In particular, we first 
want to determine the following overlap 

(J,M;Z,s = l|J,Af,A) =? (Bl) 



In order to do so, we express both states in Eq.fBTj) in terms of orbital angular momentum and spin states, which 
is pretty simple for the left hand side. The states of total angular momentum J can be written as a combination of 
states with definite orbital angular momentum I and spin s 

\J,M;l,s) = ^C(mTO^(/s)JM)|?,TO)|s,TO^) , (B2) 

where s = 1 is the spin of the vector meson, nis the z-projection of the spin, m the z-projection of the orbital 
angular momentum, M = m + iris and C is a Clebsch-Gordan coefficient. We choose the following notation for the 
Clebsch-Gordan coefficients 

{jij2,mim2\jij2,jm) = C{mim2{jij2)jm)S^,mi+-m2 ■ (B3) 

Next we turn to the state | J, M, A) in Eq. (|Bip . We want to express the helicity states of the moving system in terms of 
the spin and orbital angular momentum states. Since the spin and the orbital angular momentum are not conserved 
quantum numbers in a relativistic framework, the helicity states will be a mixture of different states. We need the 
following relations 



l^™> - I dn\9,^)Dl*i<t>,d,0), (B4) 



47r 



and the inverse of that equation, which is 



^T.\^,rn)^^DUcj,,0,O). (B5) 
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At the same time we can write for a helicity state moving along the z axis |z, A) 

|z,A) = |z)|s,A), (B6) 

where |s, A) is the usual spin state with — A. Although the spin is not a conserved quantum number, it coincides 
with the helicity state in the rest frame of the particle. Since one can not produce any orbital angular momentum 
along the direction of motion, after a boost the z projection of the total angular momentum is still given by the spin 
projection of the particle, which is the same as the helicity. Thus, we can use the above decomposition. Next we 
apply the rotation operator ?/((/), 0, 0) to the state. After the rotation the spin and helicity states will not be the same 
anymore, but the connection is given by the Wigner rotation functions. We have to rotate each factor on the right 
hand side of Eq. (IB6P separately, which gives 



|0,0,A) =[/((^,0,O)|z,A) =E|0,0)i?i„^^(0,0,O)|l,m5). (B7) 
Applying the projection operator (see [s^ or js^]) on definite total angular momentum states, we get 



|J,M,A) = \/^^ / i^i;A(0,e,O)|0,,/),A)dl] 



ms ,^ ,771 

We use the following relation for the Wigner rotation functions 



DLn{R)Di'n'iR) ^ E C (mm' ij f ) J M)) Di^{R)C inn' Uj') J N), (B9) 

J,M,N 



which yields 



J,M,A)= E ^^^^47^/<*A(</',^,O)D5:+„^,;,(0,0,O)df7|Z,m)|ml) 



ms .l,m.l' 

C{msm{ll)l'ms + m)C{OX(ll)l' X) 



E Y V "47!^ 2TrSM,m+ms J diix{x)di.j^ix)dx\l,m)\l,ms) (BIO) 



7ns ,l,m,l 

Cimsm{ll)l'ms + m)C{OX{ll)l' X) 



E y|^t^("^s(M - ms)ill)JM)CiOXill)JX)\l, M - ms)|l, ms) 

l,ms 



Therefore, we get from Eq. (jB2p and Eq. ljBlGIl 



(J, M; 1, 1| J, M, A) = y |^C(GA(n) JA) , (Bll) 

where we used 

E CimsiM -ms)ill)JM)C{{M' ~ms')m,,ill)J'M') ^ 6jj>6mm' ■ (B12) 

ms -m^/ 

Now we will connect the helicity projectors to angular momentum projectors. In order to do so we notice that 
M^j = (J, M, A|r| J, M, A) = E^-^' ^1-^' 11^1-^' ^' 1) ^) 

1,1' 

Hi] + 1 Ml]' + 1 1 _ _ (-1^-'^"^) 

= E(^> ; I', 1|T| J, M; 1) VL^-A^ic(OA(n) JA)C(OA(/'l) JA) , 
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where we used that also the orbital angular momentum states build a complete basis. By building quotients of the 
respective amplitudes, we can pin down constraints. If we only consider — l"*" and therefore only deal with s- and 
d-waves, we know for all possible combinations of I 

Mn Mn Mn Mm 



Therefore, we can use 



Mi_i M_i_i M_ii Af_io Mo_i 

M+ = Mil + Mi_i = 2Mii 

M+ = ^(Mio + M_io) = V2A/10 

M+ = ^(Moi + A/0-1) = \/2Moi . 



1 



(B14) 

(B15) 
(B16) 

(B17) 



Looking up the Clebsch-Gordan coefficients of Eq. (IBISP , we get the following relations for the respective transitions 
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(B18) 
(B19) 
(B20) 
(B21) 



Calling the transitions with definite angular momentum Dab, where a,5 G {s,d} and supressing the Lorcntz indices, 
we get 



Dss = + "y|-^io + ^ 2 ' 

^'ds = Yll ^ ^^10 l=Yoi — loo , 

\/2 V2 



_ 2 2 

-Odd = ^11 i=Yw 7=Yoi 



2Yf 



00 ■ 



(B22) 
(B23) 
(B24) 
(B25) 



In principle we can multiply each of these expressions by an arbitrary normalisation constant, which we choose to 
be one, which means we use the above expressions. If we want to express our amplitude in terms of orbital angular 
momentum 



^'-'01^01 



-"^00 -'^00 



we have to solve the following equations 
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(B26) 
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Mti\ 
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M+ 




Cds 




M+i 




\Cdd) 




\MooJ 



(B27) 



The solution to that equation is 



/4 2V2 2V2 

2 -2V2 V2 -2 

2 V2 -2V2 -2 

Vl -V2 -V2 2 J 





(M^i) 












Csd 




Moi 




Cds 


) 






\Cdd) 



(B28) 



It is interesting to note, that only for a particle at rest, the Weinberg- Tomozawa term is a pure s-wave, while for 
moving particles, factors of jj reduce the s-wave part. 
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APPENDIX C: DETAILS ON THE ai CALCULATION 



In this appendix we explicitly give the coefficients Ai form Section llVDl which incorporate the nontrivial part of 
the calculation including the oi. The explicit expressions read 

fAgVuddvs , 



— (ci(Mi )i+C2(Af2 )i) + 2^ ^^p^ c^vJ4>v{M[gval + ^{fv ~ 2gv)a^ ) 



(CI) 



and 



^2 = - ^^^!^"'f' (ciKp(Affi)i - N/2Afp(Affi)3) + c2Kp(Aff )i - v^A^p(Mf )3)) 

V^-Tq V ^Pcrmrp 

, sr^ gVudgv j , w i/^ r, n <^v\ 

0y iVii-o ^ 



(C2) 



where c^p = \/2 and ckk' = ^ 1- The first part of the coefficients Ai contains the functions Aff ^, which result from 
adding the last two diagrams of Fig. [TT] 



Aff 1 = 



1 



M. 



( V2{s -ml- Mj) \ 
-{s-Ml-Ml,) 

( ^/2(s + ml- Af2) \ 

-{s + Ml-Ml,) 
^{s-ml-MD 



(C3) 



M, 



(C4) 



s - Ml, 

where VJ is the matrix resulting from Eq. (P5|) and is given by 

iJl V1211J2 ViiioJl Vi2wJ2\ 

Y J V2111J1 V22ii^2 V2iio>/i V2210J2 {C5) 

^llOl'^l Vi201^2 ^llOO'^l Vi200<^2 
\V2101'^1 V2201<^2 V2100'^l V2200«/2/ 

The remaining part of the coefficients Ai corresponds to the diagrams Fig. [Tib . [TTF . The contain the expansion 
coefficients of the scattering amplitude Mabij and are given by 

a7 = 2Af 2Miiii + V2uj^pMpMiii2 , af ^' = 2Af|..Afi2ii + \/2cj/fi^- Af^. A'/i2i2 , 



17 = (s - ml + M^)Mnii + V2yGMpMin2 , ^* = (s - + Af|..)Afi2ii + V2ViA//j^. A/1212 
"3" = , ^ Kp«r - V2Mp{2M^pMn2i + V2MpU;^pMn22) , 

1 



PcmKK'^ 



{uJKK-a'^^' - V2Afp(2A//|..Afi221 + V2MK'i^KK'Mi222) . 



TTp 



-^^Kpaa" - 2sA/p(V2tJ^pAfii2i + MpMn22)) , 



KK* 



crairp \ 
1 



PcmKK'V^ 



-.{ujKK'a^'^' - 2sMk'{V2uokk'Mi221 + Af^-AW)) 



(C6) 
(C7) 

(C8) 
(C9) 

(CIO) 
(Cll) 
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